§ 


■ ,  r  r 

FlB  1  1  1981 


n  ) 


MAGNETIC  ORDER  AND  ENERGETfCS  OF  DIPOLAR  COUPLING  IN 
MAGNETIC  SUPERCONDUCTORS 

Martha  Red!  and  P.  W.  Anderson* 

Princeton  University** 

Princeton,  New  Jersey  08544 


ABSTRACT 

The  exact  (not  just  nearest  neighbor)  dipolar  coupling 
energy  densities  for  the  observed  structures  of  the  rare 
earth  Chevrel  compounds  are  calculated  by  the  method  of 
Luttinger  and  Tisza.  The  dipolar  coupling  energy  density 
for  the  most  probable  spin  configuration  is  comparable  to  the 
observed  magnetic  transition  temperature  but  the  most  probable 
ground  state  is  not  found  experimentally.  The  discrepancy 
between  the  observed  magnetic  ground  state  and  that  predicted 
from  dipole  coupling  may  arise  from  conduction  electron 
effects  or  possibly  some  strong  crystal  field  effect  and 
should  be  included  in  any  electronic  theory  of  the  supercon¬ 
ductive  state  below  T„. 
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I.  Introduction 

The  antiferromagnetic  interaction  energy  of  the  rare  earth  Chevrel'  compounds  is  calculated 
on  the  basis  of  simple  dipolar  coupling.’  The  compounds  are  of  special  interest  because  they  are 
also  superconductors  and  the  connection  between  their  superconductive  and  magnetic  proper¬ 
ties  is  the  object  of  much  current  experimental  and  theoretical  '*  ''  '®  effort. 

If  the  antiferromagnetic  transition  temperature  could  be  accounted  for  by  dipole  ccupling  of 
the  rare  earth  moments,  the  connection  between  superconductive  and  antiferromagnetic  order¬ 
ing  would  be  simplified.  Essentially,  different  electrons  would  give  rise  to  the  two  effects. 

In  the  rare  earth  Chevrels  the  magnetic  moments  are  very  large,  well  separated  spatially, 
and  with  a  relatively  low  density  of  conduction  electrons.  Thus  it  seemed  possible  that  the 
molecular  field  model  of  the  transition  temperature  due  only  to  dipole  coupling  could  make  a 
significant  contribution  to  the  stability  of  the  antiferromagnetic  phases  of  the  rare  earth  Che¬ 
vrels. 

The  Chevrel  compounds  are  characterized  by  sulfur  cubes  filled  in  alternate  layers  by  Mo^ 
octahedra  or  metal,  M,,  each  4/,  or  MOf,  filled  cube  separated  from  the  next  by  an  empty 
sulfur  cube;  the  sulfur  cubes  themselves  forming  a  BCC  lattice.  This  structure  permits  the 
inclusion  of  a  great  variety  of  A/,  .  When  Afj,  is  a  superconductor,  (Sn,  Pb),  the  Chevrel  com¬ 
pound  is  found  to  be  superconducting  with  raised  transition  temperature.  is  increased  more 
for  those  of  lowest  T".,,.  If  A/,  is  not  superconducting,  (Ag,  Cu),  the  Chevrel  compound 
becomes  superconducting.  If  is  a  rare  earth  having  a  large  free  ion  magnetic  moment,  the 

Chevrel  compounds  exhibit  both  magnetic  and  superconducting  ordering.  In  particular 

• 

«« 


Also  at  Bell  laboratories.  Murray  Hill.  NJ.  07974 

The  work  at  Princeton  linivcrsit)  was  supported  in  part  b>  the  National  .SL*icncc  foundation  (irani  No  DMR  7k. 
03015.  and  in  port  by  the  U.S.  Ofticc  of  Naval  Research  Grant  No.  N(X)OI4-77-C'-()71 1 


-  2  - 


TbMOfJS^,  DyMo^fi  and  HoMOf^^  are  found  to  be  superconducting  and  at  lower  temperatures 
to  exhibit  antiferromagnetic,  antiferromagnetic  and  ferromagnetic  order. 

II.  Method  of  Luttinger  and  Tisza 

The  rare  earth  atoms  lie  on  a  nearly  cubic  lattice  (rhombohedral  angle  a  =  89.4" ).  The 
exact  mean  field  dipole  coupling  interaction  energy  for  an  ordered  array  of  moments  can  be  cal¬ 
culated  by  the  group  theoretic  method  of  Luttinger  and  Tisza".  A  simple  cubic  dipole  array  is 
denoted  by  placing  a  dipole  of  definite  moment  and  direction  at  every  lattice  point.  The  lowest 
energy  ordered  arrays  are  those  which  are  invariant  under  translations  of  the  form 
T"  =  /|(2()  +  /;(27)  +  li(2k).  are  integers,  i  J ,k  are  unit  vectors  in  the  x,y,z  directions. 

To  generate  such  oo  arrays  we  have  to  specify  8  dipoles,  one  at  each  corner  of  the  unit  cube. 
The  whole  ordered  array  is  constructed  by  the  translations  T*.  The  resulting  array  may  be  con¬ 
sidered  as  a  superposition  of  eight  oo  arrays  each  of  which  consists  of  parallel  dipoles. 

Every  array  of  class  T’  can  be  specified  by  a  set  of  24  members,  e.g.,  w,  ,  i  =  1,2,  ...  24. 

Every  set  of  24  real  numbers  defines  a  T-  array  and  denotes  a  point  in  a  24  dimensional  vec¬ 
tor  space  R.  Arrays  of  constant  dipole  strength  form  a  16  dimensional  hypersurface  in  R. 

To  compute  the  energy  of  a  particular  array  M,  it  is  necessary  to  know  the  field  generated 
by  M  at  all  the  lattice  points.  The  field  /  of  a  dipole  ?  at  a  point  F  is 

/  =  [3F  •  F)  — r’]  /  r, linear  in  fi. 

The  field  generated  by  M  will  have  the  same  symmetry  as  the  array.  The  set  of  vectors 
representing  the  field  at  the  lattice  points  will  also  correspond  to  a  vector  in  the  space  R,  and 
will  be  denoted  by  F.  F  may  be  written  F  =  FM . 

Luttinger  and  Tisza  show  that  any  array  of  the  class  T’  may  be  written  as  a  linear  combina¬ 
tion  of  the  24  basic  arrays  (Fig.  I); 

I  =  1,2  ■  •  •  8. 

These  basic  arrays  form  a  complete  set  of  24  orthogonal  vectors 


^  and  >2  are  obtained  from  ^  through  cubic  rotations.  The  field  corresponding  to  M  is 

^  =  2  («-  /.  ^  +  b.f,Y^  +  c,  f,  Z^) 

1=1 

and  the  energy 

U  =  (a,-  +  b,-  +  C,-)  /, 

1=1 

where  /,  is  the  field  at  a  lattice  point  due  to  ordered  dipoles  \’-  .n  infinite  array  Z^.  The  com¬ 
putation  of  the  energy  of  any  array  is  reduced  to  the  knowledge  of  the  characteristic  values  /, . 

The  results  of  numerical  computation  of  {/, }  for  the  24  infinite  arrays  X,.Yj,  Z,  from  Table 
II  of  Luttinger  and  Tisza  in  units  of 

/,  =  0 
/,=  -9.687 
/j=  4.844 
/^=  4.844 
/,=  5.351 
/(,=  -2.676 
/,=  -2.676 
/«  =  0 

From  this  table  the  lowest  energy  for  the  simple  cubic  lattice  is  ~f  =  -2.676  The 

corresonding  array  is  05^-!-  6s21i'*-  ^5^  ~  '•  If  the  dipoles  were 

free  to  rotate  they  would  prefer  this  ordering  with  transition  temperature  r^fZs)  in  Table  I.  A 
spherical  sample  is  assumed  unless  otherwise  specified. 
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III.  Coideasatioi  Eaergy  for  Aiteroatiog  Plaaet  of  Parallel  Dipoles 

The  powder  neutron  diffraction  scans  of  Moncton,  et  al^  for  DyMo^  5$  and  ThMo^  St  below 
Tii  show  a  simple  antiferromagnetic  ordering  with  alternating  (100)  planes  of  parallel  dipoles, 
the  moments  being  alternately  parallel  and  antiparallel  to  a  unique  direction. 

We  will  assume  a  temperature  independent  interaction  which  constrains  the  rare  earth  dipole 
moments  to  lie  in  planes  of  parallel  orientation,  the  orientation  direction  being  reversed  in  adja¬ 
cent  planes.  We  will  evaluate  the  dipolar  interaction  energy  for  such  an  ’alternating  planar' 
array  having  arbitrary  dipolar  orientation  angles  0  and 

The  array  corresponding  to  these  planes  of  parallel  dipoles  is  shown  in  Fig.  2.  At  lattice 
points  1,2,  3,  4,  the  dipole  moment  has  components 

rri;,  sin  6  cos  ^ 

my  sin  0  sin  ^ 

m,  «*  cos  0 

and  at  S,  6,  7,  8  components 

mj  “  —  sin  0  cos  ^  ■  — m, 
ir^'  *  —  ain  0  »n  ♦  ■  —nty 

m,'  —  cos  0  •  — 

These  are  the  24  components  of  H. 

For  arbitrary  0,^  in  terms  of  the  basic  arrays 

a  —m^Xj  +  w,  +  iHiZt 

and 

U  ’“—Vi  [mg  /s  +  tWy* /a  +  »»!,* /  J  with 
/j  “  —9.687,  ft  “  4.844,  becomes 
V  -  -2.422  (1-3  sin*  0  cos*  dl  mV/?’. 

In  particular,  if  0  *  ^  ~  0,  the  dipoles  are  perpendicular  to  the  plane,  U  ~  4.844 
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and  is  strongly  repulsive;  if  0  —  0,^  ,  the  dipoles  are  in  the  plane,  U  —  -2.42  u^/X^,  an 

attractive  condensation  energy;  and  if  0  » tan~'v^  dipoles  are  along  [  1,  1,  1  ] 

4 

and  U  «  0. 

Computation  of  the  dipole  interaction  energy  for  alternating  planes  of  parallel  dipoles  thus 
favors  alignment  in  the  plane  and  we  would  expect  this  structure  on  the  basis  of  dipole  coupling 
if  the  dipoles  were  constrained  in  planar  arrays.  This  can  be  rewritten  such  that  U  depends 
only  on  the  angle  a  between  the  direction  of  dipole  polarization  and  the  axis  normal  to  the 
planes. 

U  ~-Vi  ft  Il-3  cos*  aj  mV/?’. 

IV.  Estimation  of  for  Planar  Arrays  of  Dipoles 

For  any  ordered  array  we  can  compute  the  energy  density.  Now  is  the 

energy  of  orientation  of  each  dipole  due  to  interactions  with  all  the  others.  Equating  this 
energy  to  the  thermal  destabilization  by  fluctuations  in  9,^  directions  ■■  2(Vi  Tu)  allows  an 
estimation  of  the  temperature  required  to  destroy  the  antiferromagnetically  ordered  state.  A 
negative  Luttinger  and  Tisza  energy  density  corresponds  to  a  magnetic  transition  temperature 

-2U‘-*^/ka 

and  a  Tu  <0  means  the  transition  is  not  stabilized  by  dipole  dipole  coupling. 

For  a  —  Ttt  in  Table  I  have  been  calculated  from 
2 

--1-^4.844 

Ti/  are  computed  transition  temperatures  for  planar  arrays  of  parallel  and  antiparallel  dipoles 
with  dipole  orientation  in  the  plane.  R  and  m^  used  for  the  estimates  were  as  follows; 

0 

The  antiferromagnetic  temperatures  have  been  calculated  with  lattice  constant  R  •  6.551  A, 
the  lattice  constant  for  PbMog  The  sulfur  superstructure  provides  space  for  intercalcation 


of  diverse  metal  ions  and  the  S-S  distance  varies  very  little  among  the  Chevrel  compounds. 

O 

The  lattice  constants  for  Mj,  —  La,  Gd,  Ho,  Er  are  6.S24,  6.4S1,  6.449,  6.445,4.  As  Z  of  the 

rare  earth  increases,  R  decreases  in  the  series  measured  above.  Thus  R  (  PbMotSg  )  will 
cause  the  estimate  of  Tm  to  be  too  low  by  about  S%.  The  lattice  constants  for  Afx  ~  Dy,  Tb 
arc  also  close  to  PbMot  5*.*’ 

The  antiferromagnetic  transition  temperatures  in  Table  I  have  been  estimated  with  ti/„ 
except  for  T^,  Dy,  Ho  where  neutron  scattering  measurements  of  thmtnm  ^cre  available. 
Neutron  scattering  measures  the  interaction  of  the  neutron  and  rare  earth  magnetic  moments 
whereas  susceptibility  and  transition  temperature  (condensation  energy)  involve  the  square  of 
the  rare  earth  moment.  Thus,  Hx,  the  moment  appropriate  for  calculating  Ty  and  measurable 
by  susceptibility  experiments,  is  related  to  the  neutron  scattering  moment  by 
Mx  “  Unr — ■  For  Tb,  Dy  and  Ho  is  less  than  10%  larger  than  Mx- 

This  calculation  would  appear  to  explain  the  nature  of  the  antiferromagnetic  transition  in 
DyMo^St  and  TbMotSf.  However,  the  neutron  data  show  that  the  dipolar  axis  is  [1,  I,  1]. 
This  structure  has  zero  dipole  interaction  energy  in  the  simple  cubic  lattice. 

Luttinger  and  Tisza  point  out  that  for  nonspherical  samples,  — 16(4s'/z— must  be 
added  to  the  energy  of  the  Z  j  array.  For  a  long,  thin  needle  M)  and  the  maximum  correction, 

— ^  —  —2.094.  Any  ferromagnetic  state  will  exhibit  domain  structure  optimizing  the  dipolar 

energy  to  this  value,  and  we  list  the  corresponding  Tf  in  Table  2  for  Ho.  Exchange  or  other 
contributions  seem  to  stabilize  the  ferromagnetic  state  of  HoMotSt,  since  the  simple  cubic 
array  is  always  nonferromagnetic,  the  energy  constant  of  the  best  nonpolarized  SC  array  being 
-2.675. 

If  we  calculate  corrections  to  dipole  interaction  energy  for  rhombohedral  angle  »  89.4*  we 
find  for  Ms  and  planar,  a  orienutions  that  Ty  changes  very  little.  For  the  [1,  1,  1] 

direction  of  alternately  polarized  planes  of  parallel  dipoles  there  is  a  small  attractive  interaction 
of  about  4%  of  that  for  the  most  favorable  structure. 
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Some  other  interaction  must  stabilize  the  (1,  1,  l]  planar  structure.  It  is  not  compatible 
with  stabilization  by  any  simple  short  range  exchange  interaction. 

This  calculation  then  suggests  several  questions  which  may  be  answered  by  the  full  micros¬ 
copic  theory  of  superconductivity  for  the  Chevrel  compounds:  Why  are  the  dipoles  not  conden¬ 
sed  into  (fls  ^  Ci  What  constrains  them  to  lie  in  planes  of  parallel  polariza¬ 

tion?  What  determines  the  direction  of  polarization  to  be  along  [l,  1,  ll  and  not  in  the  plane? 
Most  puzzling  of  all,  the  Table  shows  that  the  dipolar  energy  is  just  the  right  order  of  magni¬ 
tude  to  provide  the  molecular  field,  and  also  that  except  for  the  direction  11 1  of  magnetization, 
all  of  the  observed  configurations  at  least  are  favorable  for  dipolar  interaction,  where  the  antifer¬ 
romagnetic  one  is  quite  implausible  as  the  result  of  any  isotropic  exchange  interaction. 

As  Sherlock  Holmes  said,  ’When  confronted  with  the  impossible,  one  must  accept  the 
implausible.’  We  would  almost  prefer  to  believe  that  the  crystal  field  or  some  peculiarity  of 
domain  structure  is  misleading  us  as  to  the  actual  observed  structure.  It  is  noteworthy  that 
neutron  and  Mossbauer  measurements  have  given  puzzlingly  different  values  for  magnetic 
moment  in  some  similar  compounds. 

It  is  quite  reasonable  that  there  should  be  a  very  strong  crystal  field  effect  orienting  the 
dipoles  in  the  1 1 1  direction.  Toaocommodate  the  unfavorable  dipolar  field,  the  state  could  then 
be  a  linear  combination  of  states  oriented  in  2  or  4  different  1 1 1  directions,  or  some  yet  more 
bizarre  accommodation.  A  second  implausible  possibility  is  that  the  induced  supercurrents  in 
the  Mot  5|  groups  could  be  sufficiently  strong  to  contribute  appreciably  to  the  dipolar  energy  or 
the  anisotropy. 
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